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ABSTRACT 

We present results from Chandra X-ray and Spitzer mid-infrared observations of the interacting galaxy pair 
NGC 6872/IC 4970 in the Pavo galaxy group and show that the smaller companion galaxy IC 4970 hosts a 
highly obscured active galactic nucleus (AGN). The 0.5- lOkeV X-ray luminosity of the nucleus is variable, 
increasing by a factor 2.9 to L7 x 10"*^ ergs"' (bright state) on ^ 100 ks timescales. The X-ray spectrum of 
the bright state is heavily absorbed (A^h = 3 x 10^^ cm"^ for power law models with F = L5 - 2.0) and shows 
a clear 6.4 keV Fe Ka line with equivalent width of 144- 195 eV. Limits on the diffuse emission in IC 4970 
from Chandra X-ray data suggest that the available power from Bondi accretion of hot interstellar gas may 
be an order of magnitude too small to power the AGN. Spitzer images show that 8/xm nonstellar emission is 
concentrated in the central 1 kpc of IC 4970, consistent with high obscuration in this region. The mid-infrared 
colors of the nucleus are consistent with those expected for a highly obscured AGN. Taken together these data 
suggest that the nucleus of IC 4970 is a Seyfert 2, triggered and fueled by cold material supplied to the central 
supermassive black hole as a result of the off-axis collision of IC 4970 with the cold-gas rich spiral galaxy 
NGC 6872. 

Subject headings: galaxies: active - galaxies:clusters: general - galaxies: individual (IC 4970) - galaxies: 
interactions - infrared: galaxies - X-rays: galaxies 



1. INTRODUCTION 

Supermassive black holes are known to exist at the centers 
of most galaxies. The strong correlation of the central super- 
massive black hole mass with the mass of the host galaxy's 
stellar bulge, as represented by the central stellar velocity dis- 
persion (Ferrarese & Merritt 2000; Gebhardt et al. 2000) or 
K-band absolute magnitude (Marconi & Hunt 2003), indi- 
cates that the evolution of the central black hole and that of 
its host galaxy are inextricably linked. Most recent work re- 
lating the activity of the central supermassive black hole to 
galaxy evolution has focused on the interaction of radio-loud 
active galactic nuclei (AGN) with their environment, in or- 
der to resolve the cooling flow problem in galaxy clusters and 
bring hierarchical models of galaxy formation into agreement 
with observations of galaxies at low redshift (see, e.g., Dunn 
& Fabian 2006; Best et al. 2006; Croton et al. 2006). Hai'd- 
castle et al. (2007) separated these AGN into two classes, 
low excitation and high excitation radio galaxies, based on 
their optical emission line characteristics, and suggested that 
the observed differences between the two classes are deter- 
mined by the fuel source and accretion mode of the central 
black hole. Low excitation radio galaxies tend to be massive, 
dominant elliptical galaxies, residing in the deep dark mat- 
ter potentials near the cores of large groups and rich clusters, 
and are surrounded by an ample supply of hot gas from the 
intracluster medium falling into their large gas halos. They 
show no evidence for a standard, geometrically thin, opti- 
cally thick, luminous accretion disk. Radiatively inefficient 
Bondi accretion of hot gas accreted from their surroundings 
onto the central supermassive black hole is sufficient to sup- 
ply the total energy, mechanical as well as radiative, observed 
in these systems. In contrast, high excitation radio galaxies 
(containing the most powerful radio galaxies known) exhibit 
high-excitation, narrow-line optical emission in their spectra. 
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characteristic of emission from an accretion disk. Hardcastle 
et al. (2007) argue that AGN activity in these galaxies is pow- 
ered by radiatively efficient accretion of cold gas acquired in a 
major merger with a massive cold-gas-rich partner. The high 
radio power and jets observed in high excitation radio galax- 
ies often may result from the formation of a rapidly spinning 
nucleus from the coalescence of the merging galaxies' two 
massive black holes (Wilson & Colbert 1995; Capetti & Bal- 
maverde 2005). 

Since cold-gas accretion is associated with the galaxy's in- 
teraction and merger history, and does not depend on the depth 
of the dark matter potential or the presence of a hot intraclus- 
ter medium, galaxy collisions of lower mass galaxies with 
dusty, gas-rich partners in galaxy groups also would be ex- 
pected to trigger AGN activity and black hole growth. Gravi- 
tational and hydrodynamical forces, induced by the galaxies' 
mutual interaction, funnel dust and cold gas from the gas-rich 
partner into the interacting, lower-mass companion galaxy's 
nuclear region, where it settles into an optically thick accre- 
tion disk before subsequently accreting onto the central black 
hole. As in the evolutionary model of Churazov et al. (2005), 
accretion onto the supermassive black hole is radiatively ef- 
ficient and radio emission from the AGN may be weak or 
absent. Such galaxy collisions were likely frequent at high 
redshift (z ~ 1.3), when the fraction of blue, gas-rich star- 
forming galaxies in moderately massive galaxy groups was 
high and galaxies were rapidly evolving (Cooper et al. 2006; 
Gerke et al. 2007). Thus, while perhaps less dramatic than 
radio-loud systems, cold-gas accretion in off-axis galaxy col- 
lisions and minor mergers may play an equally important role 
in the evolution of lower mass galaxies and moderately mas- 
sive galaxy groups into what we see today. 

Observations of nearby interacting galaxies in similarly 
moderately massive groups offer a unique window into the dy- 
namical processes that triggered AGN activity in lower mass 
galaxies and influenced the coevolution of the host galaxy 
with its central supermassive black hole at this earlier epoch. 
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Fig. 1. — (top) Digitized Sky Survey Bj band image, taken witli the UK 
48incli Schmidt Telescope, of the dominant elliptical galaxy NGC 6876 and 
nearby NGC 6877 at the center of the Pavo group and the large spiral galaxy 
NGC 6872 with its tidally interacting companion galaxy IC 4970 located 
8f7 to the northwest, (bottom) Zoom-in from the image in the top panel to 
highlight the interaction of IC 4970 with the spiral galaxy NGC 6872. 



when galaxies were rapidly transforming. The Pavo group 
is such a nearby (z = 0.01338) group with a cool 0.5 keV 
intra-group medium, ^13 member galaxies and velocity dis- 
persion 425 km s"' (Machacek et al. 2005) similar to the 
properties of z ~ 1 groups observed in the DEEP2 survey 
(Gerke et al. 2005). In the top panel of Figure [T] we present 
a Bj-band image of the dominant Pavo group galaxies, show- 
ing the elliptical galaxy NGC 6876 at the Pavo group center 
and the large, gas-rich tidally-distorted S AB(rs)c spiral galaxy 
NGC 6872, located 8'J to the northwest. The smaller lentic- 
ular galaxy IC 4970, that is only 1/5 to 1/10 as massive as 
the spiral galaxy NGC 6872 (Mihos et al. 1993), is located 
1 .'12 to the north of NGC 6872's center near a break ('knee') 
in NGC 6872's northern tidal arm. The small 74kms"' line- 
of-sight velocity difference' between IC 4970 and NGC 6872 
suggests that the galaxy pair forms a spiral-dominated sub- 
group, while XMM-Newton observations of a hot, gas trail 
linking, in projection, the central elliptical galaxy NGC 6876 
and NGC 6872/IC 4970 indicate that the NGC 6782/IC 4970 
subgroup has recently passed supersonically through the Pavo 
group core (Machacek et al. 2005). Thus the Pavo group, like 
groups at high redshift, may be dynamically young and the in- 
teractions between galaxies representative of those important 
to galaxy evolution at that earlier epoch. 
IC 4970 and NGC 6872 have long been known to be a 

' All Hne-of-sight velocities in this paper are taken from the CfA Redshift 
Survey by Martimbeau & Huchra. Online versions of the catalog, supporting 
software, and docu mentation on the CfA Redshift Survey are available at 
|http://cfa-www.hMYa rd.edur huchra/zcat. 



tidally interacting galaxy pair (Vorontsov-Velyaminov 1959). 
As shown in the bottom panel of Figure [T] tidal tails extend 
outward more than 30 kpc 2.') to the east and west from the 
tip of NGC 6872's spiral arms. A stellar bridge connects the 
'knee' in NGC 6872's northern arm and tidal tail to the lentic- 
ular companion galaxy IC 4970, indicating that tidal interac- 
tions between IC 4970 and NGC 6872 are ongoing. The spiral 
galaxy NGC 6872 is gas rich. From CO J=l-0 line emission, 
Horellou & Booth (1997) infer 9.6 x IO^Mq of molecular hy- 
drogen in NGC 6872's central region. Although atomic hy- 
drogen is not observed near the center of the spiral galaxy, 
1 .41 X 10'"Mo of HI gas is distributed along the full extent of 
NGC 6872's spiral arms and associated tidal tails (Horellou & 
Booth 1997; Horellou & Koribalski 2007) and 1.3 x IO'^Mq 
of HI gas is seen along the direction of IC 4970, that may be 
associated with the companion galaxy (Horellou & Koribalski 
2007). Ha emission traces ionized gas from recent star for- 
mation along NGC 6872's outer spiral arms to the ends of the 
tidal tails and from the 'knee' north in the stellar bridge (Mi- 
hos et al. 1993). The distribution of Ha corresponds closely 
to the distribution of young star clusters observed in this sys- 
tem (Bastian et al. 2005). However, Mihos et al. (2003) 
found no Ha emission from the central regions of NGC 6872 
or associated with the companion galaxy IC 4970. 

In this paper we focus on the effects of the 
IC 4970/NGC 6872 collision on the lower mass participant 
IC 4970, based on Chandra X-ray and Spitzer mid-infrared 
observations that show IC 4970 hosts a highly-obscured ac- 
tive nucleus. Results for the two massive galaxies, NGC 6872 
and NGC 6876, and the Pavo group gas from these obser- 
vations will be reported in subsequent papers (Machacek 
et al. 2007a,b in preparation). Our discussion is organized 
as follows: In ^we briefly review the Chandra and Spitzer 
observations and our data reduction and processing proce- 
dures. In ^we present our main analysis results. We discuss 
the X-ray time variability of the nuclear point source in ^3.11 
its X-ray spectral properties in 33.21 and mid-infrared flux 
densities and colors in 33.31 and show that these properties 
are consistent with the classification of IC 4970's nucleus as 
a highly obscured Seyfert 2. In 33.41 we briefly discuss the 
properties of a nearby ultra-luminous X-ray source and in 
33.5l place limits on the amount of hot gas outside the central 
region of the galaxy. In ^ we determine the mass of the 
central black hole in IC 4970 and discuss possible accretion 
modes, concluding that nuclear activity in IC 4970 is most 
likely triggered and fueled by cold gas driven into the nucleus 
during IC 4970's ongoing off-axis collision with the dust- and 
gas-rich spiral galaxy NGC 6872. We summarize our results 
in ^ Unless otherwise indicated, we quote 90% confidence 
levels for the uncertainties in spectral parameters and Icr 
uncertainties on counts and count rates. Coordinates are 
J2000.0. Adopting the flat ACDM cosmology from the three 
year WMAP results (Hq = 73 kms"' Mpc"' , n,„ = 0.238, 
Spergel et al. 2007), the luminosity distance to the Pavo 
Group (z = 0.01338) is 55.5 Mpc and 1." corresponds to a 
distance scale of 0.262 kpc. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Chandra Observations 

Our X-ray data consist of two observations of the Pavo 
galaxy group, including IC 4970, taken with the ACIS-I de- 
tector (Garmire et al. 1992; Bautz et al. 1998) on boai'd the 
Chandra X-ray Observatory on 2005 December 14-15 (OB- 
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SID 7248) and 2005 December 16-17 (OBSID 7059) with to- 
tal exposures of 35.3 and 40.7ks, respectively. The nominal 
ACIS-I aimpoint was adjusted -2.'0 in the negative Y-axis di- 
rection away from the nucleus of the spiral galaxy NGC 6872 
to prevent chip gaps from falling on the tidal features and 
diffuse emission regions of interest. The data were analyzed 
using standard X-ray analysis packages CIA03.3, FTOOLS, 
and XSPEC 11.3.0. The data were filtered to remove events 
with bad grades (1,5, 7), and those that fell on hot pixels. 
Events, flagged in VR4INT mode as having excessive flux 
in the border pixels surrounding event islands were also re- 
moved, resulting in a factor > 2 improvement in the rejec- 
tion of particle backgrounds at energies below ~ 1 keV. We 
then reprocessed the data and created response files using 
the most recent gain tables and instrumental corrections, in- 
cluding corrections for the charge transfer inefficiency of the 
ACIS-I forward illuminated CCDs, the time-dependent build- 
up of contaminants on the telescope optical filter (Plucinsky 
et al. 2003), and the secular drift of the average pulse height 
amplitude for photons of fixed energy (tgain).^ Periods of 
anomalously high and low count rates were removed from 
each observation using the script lc_clean acting on events 
in the 0.3 - 12 keV energy band from the 10 CCD, where there 
were no bright sources. This resulted in useful exposures of 
31,878s for OBSID 7248 and 38 , 364 s for OBSID 7059. 

Background files for image analysis were gen- 
erated using the 1.5 x 10^ s source free data set 
acisi_D_01236_bg_evt_010205.fits, appropriate for the 
date and instrument configuration of our observation.^ The 
background files were normalized by comparing rates in 
the 9.0-1 1.5 keV energy band, where particle background 
dominates. We found that the 9.0-11.5keV event rates 
in our observations were 18.2% (18.9%) higher in OBSID 
7248 (7059), respectively, than that found in the source free 
datasets. We thus applied an additional scaling factor of 
1.182 (1.189) to the source free data sets for OBSID 7248 
(7059), respectively, to bring the 9.0- 11.5keV source free 
event rates into agreement with the observations. 

2.2. Spitzer Observations 

We observed NGC 6872/IC 4970 in the mid-infrared (3.6, 
4.5, 5.8, and 8.0/im wavebands) on 2005 September 18 us- 
ing the Infrared Array Camera (IRAC; Fazio et al. 2004; 
Hora et al. 2004) on board the Spitzer Space Telescope 
(Werner et al. 2004), as part of a 432 s Spitzer IRAC obser- 
vation (AORID 14699264) of the Pavo Galaxy Group. Our 
IRAC observation of the Pavo group used two fixed positons 
to ensure coverage of both the central group elliptical galaxy 
NGC 6876 and the NGC 6872/IC 4970 interacting galaxy 
pair. We used the 12 s HDR mode and a 9-position, large- 
scale random dithering pattern to minimize latency artifacts, 
while providing sufficient redundancy for reliable removal of 
cosmic rays and scattered light artifacts. This resulted in an 
effective integration time at the position of IC 4970 of 108 s in 
the 3.6 and 5.8/im wavebands and 216 s in the 4.5 and 8.0/im 
wavebands. We used Basic Calibrated Data reprocessed by 
the Spitzer Science Center S14.0 pipeline on 31 May 2006 to 
include the super-boresight correction for improved astrome- 
try. The Basic Calibrated Data frames in all four IRAC wave- 
bands were cleaned using the script cosmetic.pl in the Spitzer 
Science Center data analysis software suite MOPEX (version 

^ See |http://cxc. harvard.edu/contrib/alexey/tgai n/tgain.html I 
' see[http://cxc.harvard.edu/contrib/maxim/acisbg] 
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Fig. 2. — (top) Spitzer IRAC 3.6^im closeup image of IC 4970, tracing the 
distribution of old (MO) stars. Contour levels are 0.02 and 0.12 mjy arcsec"^. 
The image field of view is 32"4 X 5 1 "4 with 1 pixel = 0" 86 X 0" 86. (bottom) 
Chandra 0.5 — 2 keV image of IC 4970 with IRAC 3.6/im contours from the 
top panel superposed. The X-ray emission is concentrated in two X-ray point 
sources, one at the nucleus and one located 2"6 (0.68 kpc) to the north of the 
nucleus. 1 pixel = 0"492 X 0"492. 



030106) to mitigate muxiplexer bleed, and to correct for col- 
umn pulldown and dropouts. We then used the 'post-Basic 
Calibrated Data' processing suite IRACproc 4. 1 .2 (Schuster 
et al. 2006) with linear interpolation and 0"86 pixel size to re- 
ject outliers and prepare mosaics, registered on the same sky 
grid, in each waveband. A small, residual background gradi- 
ent persisted across the mosaics in each channel. We modeled 
this background gradient using the SExtractor 2.5.0 software 
package (Berlin & Arnouts 1996). All of the resulting back- 
ground maps were inspected and found free of anomalies that 
would contaminate the photometry. We then subtracted the 
background map for each channel from the data to create the 
final background-subtracted mosaics. 

3. DATA ANALYSIS 

In the top panel of Figure|2] we show a Spitzer IRAC 3 .6/im 
image of IC 4970, tracing the stellar distribution in the galaxy. 
For reference, 0.02 and 0.12 mJy arcsec"^ contour levels for 
the 3. 6/im emission are also shown. The 3. 6/im emission seen 
in the lower southeast corner of the image is from the north- 
ern spiral arm of NGC 6872. In the bottom panel of Fig. |2]we 
show the 0.5-2keV Chandra X-ray image of the same field, 
using the superposed 3.6//m contours from the top panel to 
delineate the stellar extent of IC 4970. X-ray emission from 
IC 4970 is concentrated in two point sources within the cen- 
tral 1 kpc (4") of the galaxy with one point source located at 
the nucleus (RA = 20''16"'57''.56, DEC = -70°44'59."9), and 
one 0.68 kpc (2"6) to the north (RA = 20''16'"57^6, DEC = 
-70°44'57"3). The analysis regions we use to study these 
sources and diffuse emission in IC 4970 are listed in Table[T] 
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Fig. 3. — Chandra X-ray light curves for an elliptical region with semiaxes of 2"2 and 1"7 centered on the nucleus of IC 4970 in 3.2 ks bins, showing a factor 
of three increase in mean count rate over the ~ 136 ks time period between Obsid 7248 (left panel) and Obsid 7059 (right panel). The mean time zero in both 
panels is measured from the start of Obsid 7248. 



3.1. Nuclear X-ray Variability 

The X-ray luminosity of IC 4970 is dominated by emis- 
sion from the nuclear point source. Enhanced nuclear activity 
is often found as a result of gas rich mergers, either in the 
form of a circumnuclear starburst or increased AGN activ- 
ity. To explore the nature of the nuclear source in IC 4970, 
we constructed light curves in the 0.5- lOkeV energy band 
for our two observations, that were taken 136ks apart, using 
an elliptical region with semiaxes of 2. "2 and 1"7 (denoted 
'nucleus' in Table [Til and 3.2 ks binning to reduce statisti- 
cal scatter. Errors were calculated using the Gehrels method 
for low count rate data (Gehrels 1986). Light curves for 
the nucleus, calculated with and without background sub- 
traction, showed no significant difference. Figure [3] shows 
that the nucleus is highly variable, with the mean count 
rates, (6.5 ± 1.8) x lO-^cntss"' for Obsid 7248 (left panel) 
and (19.2 ±2.5) x lO'^cntss"' for Obsid 7059 (right panel), 
increasing by a factor ~ 3 over the 136ks interval between 
the observations. To check for any systematic variation in 
count rate between the two observations, we used the same 
energy band and time binning to compute the light curve for a 
30" circular region outside IC 4970 (denoted TGM' in Table 
[TJ. We find no statistically significant time variation between 
the two observations for region TGM' . Scaling the observed 
IGM count rate by the ratio of the solid angles for the two 
regions ('nucleus' and 'IGM'), we find that the mean IGM 
background count rates forlC 4970's nucleus, (1.86 ±0.62) x 
10"^cntss"' and (2.2 ±0.75) x 10"^cnts s"' for the two obser- 
vations, agree within their Icr uncertainties. 

Short time scale variability is not expected for circumnu- 
clear starbursts. Thus the observed variability strongly sug- 
gests that the nucleus of IC 4970 is an AGN. 

3.2. Nuclear X-ray Spectra and Luminosities 

We extracted the spectrum of IC 4970's nucleus in the same 
elliptical region used to construct the light curves, and used a 
local background annulus concentric with the nucleus with 
inner and outer radii of 4."9 and 19."7, respectively, (denoted 
'nucleus bkg' in TablefTJ to subtract the background. 

3.2.1. Bright State 

We found 683 ±26 net source counts in the 0.5- lOkeV en- 
ergy band for the bright state (Obsid 7059). Bright unresolved 



sources observed with ACIS may be subject to pileup, where 
two or more photons incident on the same detection cell are 
recorded as a single event with a pulse height amplitude that 
is roughly the sum of the individual event pulse height ampli- 
tudes (Davis 2001). From the observed mean count rate per 
3.2s frame (0.06 cnts frame"') for the nucleus in the bright 
state, we used PIMMS to estimate a pileup fraction of ~ 3% 
in our data, such that pileup should not significantly affect our 
results. We grouped the spectrum to ensure a minimum of 30 
counts per bin. As shown in the left panel of Figure |4] the 
spectrum is highly absorbed and shows a clear 6.4 keV Fe Ka 
line. 

We model the spectrum over the 3- lOkeV range using a 
power law model with an added narrow Gaussian to model 
the Fe Ka line. The model is corrected for absorption us- 
ing photoelectric cross sections from Morrison & McCam- 
mon (1983). For the narrow Gaussian line profile, we fix the 
central line energy at 6.4 keV, and redshift at z = 0.015728. If 
we allow all other parameters to vary freely, we find a best fit 
hydrogen absorbing column of 2.9;!;[ ^ x lO^-' cm~^ . Although 
the central value for the photon index from this fit (F = 1 .6) 
is consistent with that expected for an AGN, F is poorly con- 
strained. We thus also fix the photon index F in our mod- 
els to lie in the range 1.5-2.0 typical for highly obscured 
AGN (Cappi et al. 2006). The results of our spectral fits are 
summarized in Table |2] For F = 1.5-2.0, we find column 
densities of 2.8 - 3 x lO^-' cm"^ with 10% uncertainties and Fe 
Ka line equivalent widths of 144- 195 eV. The intrinsic X- 
ray luminosity of the nucleus in the bright state is high, with a 
0.5-lOkeV luminosity ofl.l'^^l^ x 10"*^ erg s"' , where the er- 
rors reflect the uncertainty in the photon index in the spectral 
models. 

3.2.2. Faint State 

The spectrum of the faint state (Obsid 7248) is shown in the 
right panel of Figure |4] There are 196 ± 14 net source counts 
for the faint state in the 0.5- lOkeV energy band with only 
15.6 ±3.9 net source counts below 2keV. We grouped the 
spectrum using a predefined grouping for very faint sources 
that results in channels of approximately constant logarith- 
mic width, and restricted our spectral fits to the 2- lOkeV 
energy band. We again model the spectrum using absorbed 
power law models with fixed photon index F in the range 
1.5-2.0. The x^/dof in the model fits to the faint state fa- 
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Fig. 4. — (left) Spectrum of the nucleus of IC 4970 during Obsid 7059 (bright state). The spectral model (solid line) is an absorbed power law plus Gaussian, 
with hydrogen column 3 X 10^' cm"^ , photon index 1.7, and a narrow Gaussian with line energy and redshift of 6.4 keV and 0.015728, respectively, (right) 
Spectrum of the nucleus of IC 4970 during Obsid 7248 (faint state). The spectral model (solid line) is an absorbed power law with photon index 1 .7 and hydrogen 
column 3.1 X 10^^ cm"^ . Please note the difference in vertical scales, and each major tick mark in the right residual panel is 2 X 10"^. 




vor Steeper T ~ 1.7-2.0, as expected for a Seyfert 2 nu- 
cleus. The hydrogen absorption column for these models is 
(2.9-3.4) X 10^^ cm"^ , with ±20% uncertainties from the 
spectral fits, in agreement with the spectral model fits for the 
bright state (See Table |2]l. The 0.5- lOkeV X-ray luminos- 
ity of the nucleus in the faint state is 0.6;!;j]'j x lO^'^ergs"' , 
where the uncertainties are due to the uncertainty in the spec- 
tral model r. 

Thus, with these data, we see no significant difference be- 
tween the bright and faint states in either the power law slope 
or hydrogen absorption column. We find only that the X-ray 
luminosity is reduced by a factor 2.9. 

3.2.3. Soft Emission in the Nucleus 

There are a total 41 ±6.4 source counts from the nuclear 
region in the 0.5-2keV energy band for the combined ob- 
servations. If we assume that this soft emission is solely due 
to thermal emission from hot gas in the nuclear region out- 
side the obscuring torus, we can place an upper limit on the 
amount of hot gas within ^ 0.5 kpc of the nucleus. Using 
Anders & Grevasse (1989) abundances and assuming solar 
metallicity, the spectrum of this soft emission is consistent 
with an absorbed thermal APEC model (Smith et al. 2001) 
with absorbing column A^h = (8 ±4) x 10^' cm~^ and temper- 
ature kT = 0.6!j]-^keV (x^/dof = 6.1/7), and an estimated 
0.5- lOkeV intrinsic luminosity of 9 x lO-'^ergs"' . From 
the spectral model normalization and assuming uniform fill- 
ing, we estimate an electron density and hot gas mass in the 
nuclear region of < 0.2 cm"'' and < 3 x 10''Mo, respectively. 

3.3. IRAC Mid-Infrared Colors of the Nucleus 

Heavily absorbed nuclear regions of galaxies are typically 
bright in the mid-infrared. Figure |5] shows a composite IRAC 
mosaic of the NGC 6872/IC 4970 system. Blue and green 
denote 3.6/im and 4.5/xm emission, respectively, that trace 
old stars in the galaxies. The 8/xm emission is a combi- 
nation of starlight and 'nonstellar' emission sources such as 
warm dust and AGN. We model the stellar contribution by 
averaging background-subtracted 3.6/im and 4.5/im emission 
maps that have each been aperture corrected and scaled to the 
8/im waveband using the mid-infrared colors of MOIII stars. 
We then subtract this model from the background-subtracted 
8/xm mosaic to produce a map of 'nonstellar' 8/im emission, 
shown in red (Pahre et al. 2004; Machacek et al. 2007a, in 



Fig . 5 . — Spitzer mid-infrared composite image of NGC 6872 and IC 4970 
taken in the IRAC wavebands. Blue (green) shows the 3.6(4.5)/.im emission 
from old stars, respectively, and red shows the 8/xm nonstellar emission (see 
Machacek et al. 2007a, in prep.). North is up and east is to the left. 



preparation). Figure |5] shows that the nonstellar 8/zm emis- 
sion in IC 4970 is concentrated in the inner 4" -5" (~ 1 kpc) 
of the galaxy, with the peak of the nonstellar emission at 
the nucleus of IC 4970, qualitatively consistent with a heav- 
ily obscured active nucleus (Pahre et al. 2004). To check 
whether the 8/im nonstellar emission is extended or consis- 
tent with an unresolved point source, we extracted 2 "2 wide 
vertical and horizontal cuts through the 8/im nonstellar sur- 
face brightness distribution at the nucleus. The full width at 
half-maximum of the resulting projected surface brightness 
distribution in each case was ^ 3.", a factor 1.5 larger than 
that of the IRAC point spread function in this waveband, sug- 
gesting part of the emission arises from a roughly spherical 
distribution of warm dust in the central region of the galaxy 
outside the luminous accretion disk of the AGN. The bright 
patches of nonstellar 8/^m emission from polycyclic aromatic 
hydrocarbons (warm dust) along NGC 6872's northern arm 
and at the 'knee' and 'bridge', that trace the collision path of 
IC 4970 through NGC 6872, also show that IC 4970 has re- 
cently tidally interacted with a region in the spiral that is rich 
in cold gas and dust, such that the spiral galaxy NGC 6872 
could be one source for the obscuring material now found in 
the less massive, tidal perturber IC 4970. 
We perform fixed aperture photometry on the background- 
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subtracted mosaiced-images in all fourlRAC wavebands (3.6, 
4.5, 5.8, and 8.0/im), using a 2"44 circular aperture and back- 
ground annulus with inner, outer radii of 2 "44, 7. "32, respec- 
tively, centered on the nucleus of IC 4970. The data were 
aperture corrected using point source aperture corrections for 
each IRAC waveband appropriate for the chosen aperture and 
background geometries, and the resulting flux densities con- 
verted to Vega relative magnitudes. For completeness, the 
aperture correction factors and Vega relative zeropoints used 
for each IRAC channel are listed in Table [3] Our results for 
the IRAC -band flux densities (magnitudes) of IC 4970's nu- 
cleus are given in Table H] 

We find IRAC mid-infrared colors of [3.6]- [4.5] = 0.17 
and [5.8] -[8.0] = 1.59. These mid-infrared colors fall out- 
side the mid-infrared color selection criteria used by Stern 
et al. (2005) to separate active galaxies from normal galax- 
ies. The [3.6] - [4.5] color of IC 4970's nucleus is too blue, 
falling below the color-color selection 'wedge' shown in their 
Figure 1, but the [5.8] - [8.0] color of IC 4970's nucleus is 
far redder than the zero redshift endpoint of the evolution- 
ary track for the normal SO/Sa galaxy NGC 4429, also shown 
in that figure. However, the Stern et al. selection 'wedge' 
was based on a photometrically identified sample of bright 
quasars and Seyfert 1 galaxies, rather than highly obscured 
Seyfert 2s. Given that 60% of the narrow line AGN iden- 
tified in their sample also fall outside the color-color selec- 
tion 'wedge', these criteria likely do not apply to IC 4970. 
Lacy et al. (2004) used quasars from the Sloan Digital Sky 
Survey to empirically define a selection region in the mid- 
infrared color-color space defined by the 8.0/im/4.5/^m and 
5.8/im/3.6/im flux density ratios. Objects red in both colors 
are likely AGN. These mid-infrared colors for the nucleus of 
IC 4970 ( log(55.8/53.6) = 0.13 and log(S&,Q/S4.5) = 0.6) are 
both red and are consistent with the colors of obscured AGN 
candidates in the Lacy et al. sample (See Fig. 1 in Lacy 
et al. 2004). 

Thus both the X-ray properties (spectrum and variability) 
and mid-infrared colors suggest IC 4970 hosts a highly ob- 
scured Seyfert 2 nucleus. 

3.4. Near-Nuclear X-ray Source 

We use a l."3 circular source region and the same annu- 
lar background region that was used for the nucleus ('north 
source' and 'nucleus bkg', respectively, in Table[T]i to measure 
the net X-ray source counts in the combined 70.2 ks observa- 
tion from the near-nuclear X-ray source, located 0.68 kpc to 
the north of the AGN (see the right panel of Figure|2]i. Since 
the near-nuclear X-ray source is separated from the nucleus by 
only ^ 2."6 and the nucleus is located 3.'19 from the ACIS- 
I aim point for this observation, the AGN and near-nuclear 
point source are resolved in the 0.5-2keV energy band, but 
are not resolved in the 2- lOkeV energy band. We find a total 
of 37.5 ±6.2 net source counts in the 0.5-2keV energy band 
for the northern, near-nuclear X-ray source and, to avoid con- 
tamination from the nucleus, use only the 0.5-2keV count 
rate to determine its luminosity. Assuming Galactic absorp- 
tion of 5 X 10^° cm"^ , and absorbed power law spectral mod- 
els with photon indices F = 1.7 ±0.3, typical for luminous 
X-ray point sources (see, e.g. Colbert et al. 2004; Swartz 
et al. 2004), we find a 0.5- lOkeV luminosity for the near- 
nuclear X-ray point source of 3.5;!;[ f x lO-'^ erg s"' , where the 
errors reflect both the range of photon indices considered in 
our spectral models and the la statistical uncertainty in the 
0.5-2keV count rate. Thus the near-nuclear X-ray point 



source is an ultra-luminous X-ray source. Our estimate for the 
luminosity of the near-nuclear X-ray source should be consid- 
ered a lower bound on its intrinsic luminosity, since 8/im non- 
stellar emission in the center of IC 4970 suggests that obscur- 
ing material may extend outside the nucleus (see ^3. 3b . and 
ultra-luminous X-ray sources, on average, tend to be moder- 
ately absorbed (A^h 2 x 10^' cm"^ ; Swartz et al. 2004). 

3.5. Diffuse Emission Outside the Central Region 

To estimate the amount of hot gas in IC 4970 outside the 
nucleus, we measure the X-ray emission from an elliptical 
region, based on the distribution of starlight in the image of 
IC 4970 shown in the right panel of Figure [T] excluding a 
central ellipse containing the two point sources near the nu- 
cleus (denoted 'full extent' and 'central' in Table [TJ. The lo- 
cal background is again the nearby 30." circular region 'IGM' 
in Table [U X-ray emission outside the central region is weak 
with only 28.3 ±6.4 net source counts in the 0.5-2keV en- 
ergy band, for the combined 70.2 ks observation. We estimate 
the amount of this emission that is the result of low mass X- 
ray binaries in IC 4970 using the low mass X-ray binary cu- 
mulative luminosity function for early type galaxies from Gil- 
fanov (2004), normalized by the stellar mass determined from 
K-band photometry of the 2 Micron All Sky Survey image of 
the same region. We model the excess 0.5-2keV emission 
with a Raymond-Smith thermal plasma model, assuming a 
gas temperature of 0.54 keV and metallicities of 0.4- IZq, 
typical for the interstellar medium of early type galaxies (see, 
e.g. Machacek et al. 2004, 2006), and Galactic absorption 
of 5 X 102"cm-2 (Dickey & Lockman 1990). We find a total 
0.5-2keV luminosity outside the central region of IC 4970 
of 10^*^ erg s"' , with 6.8 x 10^** erg s"' (66%) from the thermal 
component and 3.5 x 10""^ erg s"' (34%) from low mass X-ray 
binaries. Using the Raymond-Smith spectral model normal- 
ization to estimate the mean electron density and gas mass 
outside IC 4970's central region, we find (for uniform fill- 
ing) mean electron densities of (2 - 3) x lO"-' cm"-' and hot gas 
masses of '--^ 1 - 2 x lO^M©, where the ranges reflect uncer- 
tainties in the assumed metal abundance and geometry. This 
is a factor > 50 less than the amount of cold HI gas that may 
be associated with the galaxy (Horellou & Koribalski 2007). 
Thus either the gravitational potential of IC 4970 is too shal- 
low to hold a significant hot gas halo or, alternatively, the hot 
gas has been stripped during the collision of IC 4970 with 
the spiral galaxy NGC 6872 and/or by its interaction with the 
Pavo group intra-group medium during the supersonic pas- 
sage of the NGC 6872/IC 4970 galaxy pair through the Pavo 
group core. 

4. FUELING THE CENTRAL BLACK HOLE 

4. 1 . Determining the Black Hole Mass 

The mass of the central black hole in galaxies is tightly 
correlated with the absolute K-band luminosity of the galaxy 
bulge, 

logMBH = 8.21 + 1.131og(L,f- 10.9) (1) 

where Mbh is the mass of the black hole in solar masses 
and and Lk is the absolute K-band luminosity of the galaxy's 
bulge in solar luminosities (Marconi & Hunt 2003). For early 
type galaxies, like IC 4970, the bulge luminosity in eq. [T]is 
replaced by the total K-band luminosity of the galaxy as a 
whole. Thus we use photometry of the 'full extent' source 
region from Table [T] with the 'IGM' region as a local back- 
ground, on 2.2/im images from the 2 Micron All Sky Sur- 
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vey to measure the K-band luminosity of IC 4970. Since it 
is likely that IC 4970 and NGC 6872 lie near the core of the 
Pavo group (Machacek et al. 2005), such that the luminos- 
ity distance is 55.5 Mpc, we find logiLx/LQ) = 10.46 and the 
central black hole mass of 5 x IO^Mq. If instead we use the 
B-band mass-to-light ratio (Mdyn/ie = 10 ) for IC 4970, that 
was needed in simulations to reproduce the tidal features ob- 
served in the spiral galaxy NGC 6872 (Mihos 1993; Horellou 
& Koribalski 2007), we find a dynamical mass of ^ 10"Mq 
within a B-band effective radius of 11" (3kpc), implying a 
line-of-sight velocity dispersion of 131 km s"' . Then from 
the Mbh-o" relationship (Gebhardt et al. 2000; Ferrarese & 
Merritt 2000), we infer a central black hole mass for IC 4970 
of 2 - 3 X IO^Mq, a factor ~ 2 smaller than that obtained from 
the Lk-Mqh relationship, but consistent within its ±0.5 dex 
scatter (Hardcastle et al. 2007). 

4.2. Mass Accretion Rates and Modes 

From the X-ray properties of the nuclear region of IC 4970, 
we can estimate the accretion power available from hot gas in 
the nuclear region, and compare the available accretion power 
from hot gas to the luminosity of the nucleus to constrain the 
accretion mode of the AGN. For a 5 x 1O''M0 black hole the 
Eddington luminosity is Lndd = 6.3 x lO'^^ergs"' , such that 
for the bright state, Lx/L-em = 2.7;!^ j x 10"^, consistent with 
Lx I i'Edd observed for heavily absorbed, accretion powered ra- 
dio galaxies (Evans et al. 2006). Thus, even in the bright 
state, the radiative power is sub-Eddington. The mass accre- 
tion rate needed to produce the observed X-ray luminosity is 
given by ffi = Lx/{r]c^), where 77 is the radiative efficiency and 
c is the speed of light. Assuming a 10% radiative efficiency, 
the mass accretion rate needed to power the AGN in IC 4970 
is rii = 1 .9^0 3 X 10^^ g s"' . This should be considered a lower 
limit on the required mass accretion rate, since for accretion 
powered AGN, the bolometric luminosity is likely a factor 
^3-10 higher (Elvis et al. 1994) than Lx, and radiative effi- 
ciencies < 0.01 may be more common (Evans et al. 2006). 

The accretion rate available from the Bondi accretion of hot 
gas onto the central black hole is given by 

MBondi = 3^ (2) 

where G is the gravitational constant, Mbh is the black hole 
mass, Cj is the sound speed and p is the gas density at the 
Bondi radius ka = IGMqu/c^ (see, e.g., Allen et al. 2006). 
In our case the Bondi radius is small (^ 3pc), so we use the 
limits on the hot gas near the nucleus derived from assuming 
all of the excess 0.5-2 keV emission in the nuclear region was 
thermal (see ^3.2l i. For hot gas with temperature kT keV 
and electron density n^ ^ 0.2 cm"-' accreting onto a 5 x lO^M© 
black hole, we find a Bondi accretion rate of 8.4 x 10^" g s"'. 
This is less than 5% of the lower limit on the mass accretion 
rate ni needed to power the AGN. 

The available Bondi accretion power is very sensitive to the 
temperature of the hot gas through the dependence of eq. |2] 
on the gas sound speed, with lower gas temperatures allow- 
ing higher accretion rates. However, even for the 90% confi- 
dence lower limit on the gas temperature {kT = 0.2 keV) for 
hot gas in the nuclear region, the available mass accretion rate 
(6.8 x 10^' g s~') is still a factor of three too small to power the 
observed X-ray luminosity and likely an order of magnitude 
too small to provide the total bolometric power for the AGN. 
Furthermore, if we have overestimated the black hole mass 



by a factor two, as suggested by previous numerical simula- 
tions of the interaction of IC 4970 with the larger spiral galaxy 
NGC 6872 (Mihos et al. 1993, Horellou & Koribalski 2007), 
the Bondi accretion rates quoted above would be reduced by 
an additional factor of four. 

Thus the most likely explanation for the nuclear activity 
in IC 4970 is that the active nucleus was triggered by the 
accretion of cold (HI) gas and dust driven into its center 
as a result of IC 4970's ongoing interaction with the spi- 
ral galaxy NGC 6872. Since numerical simulations of the 
NGC 6872/IC 4970 system by Horellou & Koribalski (2007) 
suggest that IC 4970 will accrete significant amounts of mat- 
ter during this interaction, some, if not all, of the dust and 
cold gas fueling IC 4970's nucleus may have been acquired 
by IC 4970 during its passage through the northern spiral 
arm of its much larger, gas-rich spiral partner. Such interac- 
tions would have been frequent between galaxies in groups at 
high redshift, when the fraction of gas-rich spirals was high, 
thereby promoting AGN activity and black hole growth in 
lower mass galaxies, that are unable to retain significant hot 
gas halos because of their shallow gravitational potentials. 

5. CONCLUSIONS 

In this paper we use Chandra X-ray and Spitzer IRAC mid- 
infrared observations to probe nuclear activity in IC 4970, 
the interacting companion galaxy to the large spiral galaxy 
NGC 6872 in the Pavo galaxy group. We find the following: 

• X-ray emission from the nucleus of IC 4970 is time 
variable with the mean 0.5- lOkeV count rate increas- 
ing by a factor of three on ^ 100 ks time scales. This 
strongly suggests IC 4970's nucleus hosts an AGN 
rather than a compact nuclear starburst. 

• The X-ray spectrum of IC 4970's nucleus in the bright 
state shows a clear Fe Ka line and is well described 
by an absorbed power law plus narrow Gaussian line 
model with absorption column A^h = 3 x 10^^ cm"^, 
fixed photon indices F = 1.5-2, and an Fe Ka line 
equivalent width of 144- 195eV. The 0.5- lOkeV X- 
ray luminosity of the nucleus in the bright state is 
1.7+03 ^ 10**^ ergs"' and in the faint state is 0.6:l:oi x 
lO'*^ erg s"' . With these data, we find no significant dif- 
ference in the power law slope or intrinsic absorption 
between the bright and faint nuclear states. 

• Nonstellar 8/im emission is concentrated in the cen- 
tral - 1 kpc (4" -5.") of IC 4970, consistent with high 
obscuration in this region. The mid-infrared colors of 
the nucleus, log(55.8/53.6) = 0.13, log(58.o/'S4.5) = 0.6, 
are consistent with the mid-infrared colors expected 
for highly obscured AGN candidate sources. Thus the 
X-ray and mid-infrared properties of the nucleus of 
IC 4970 suggest that it is a highly obscured Seyfert 2. 

• Assuming that the 0.5-2keV X-ray emission observed 
in the nuclear region is all thermal, the spectrum is 
consistent with an absorbed 0.6^o4keV APEC model 
with solar metallicity, a hydrogen absorbing column of 

8 ±4 X 10^' cm"^ , and intrinsic 0.5-10 luminosity of 

9 X 10^^ erg s"' . For uniform filling, the inferred mean 
electron density and hot gas mass in the nuclear region 
are < 0.2 cm"-' and < 3 x IQ^Mq, respectively. 
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• An ultra-luminous X-ray source with 0.5- lOkeV in- 
trinsic luminosity of > 3.5;!;j 2 x 10^^ ergs"' is found 
0.68 kpc from IC 4970's nucleus. 

• Little diffuse X-ray emission is observed outside the 
nuclear region. Assuming 0.54 keV interstellar galaxy 
gas, the X-ray luminosity of diffuse gas outside the cen- 
tral 1 kpc of the galaxy is 7 x 10"'**ergs~' and gas 
mass is ^ 1 - 2 X IO^Mq. Either the gravitational po- 
tential of IC 4970 is too shallow to retain a significant 
hot gas halo, or that halo has been stripped by the inter- 
action of IC 4970 with NGC 6872 and the Pavo group 
intra-group medium. 

• From the correlation between black hole mass and K- 
band luminosity, we expect IC 4970 to host a 5 x 
1O^M0 black hole. Bondi accretion of hot gas in the 
nuclear region of IC 4970 onto the black hole can ac- 
count for only ^ 5% of the observed X-ray luminos- 
ity of the AGN, such that the dominant power source 
for the active nucleus is likely the accretion of dust 
and cold gas driven into IC 4970's nuclear region by 
the galaxy's ongoing interaction with the gas-rich spi- 
ral galaxy NGC 6872. 
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TABLE 1 
IC 4970 Analysis Regions 



Source 


Geometry 


Center 


Radii 


PA 






RA, DEC 


ai'csec,arcsec 


deg 


Nucleus 


ellipse 


20 16 57.6,-70 44 59.9 


2.2,1.7 





Nucleus Bkg 


annulus 


20 16 57.6,-70 44 59.9 


19.7,4.9 




North Source 


circle 


20 16 57.6,-70 44 57.2 


1.3 




Central 


ellipse 


20 16 57.8,-70 44 59.5 


4.4,3.9 





Full Extent 


ellipse 


20 16 57.2,-70 44 59.0 


18,7.5 


96 


IGM 


circle 


20 16 48.2,-70 44 32.1 


30 




Note. — 


Columns are 


(1) region label, (2) region 


geometry, (3) re 


gion 



center in J2000.0 coordinates, (4) radius (circle), outer, inner radii (annulus), 
and major, minor semiaxes (ellipse), and (5) position angle measured counter- 
clockwise from west for elliptical regions. 



TABLE 2 

IC 4970 Nuclear X-ray Spectral Fits 



State r Nh FeEq Width x^/dof 

10-3 cm-2 eV 



Bright 



Faint 



1.5 


2.8±0.3 


144 


17.14/18 


1.7 


3 0+*-3 


171 


17.5/18 


1.8 


^•^-0.4 


179 


I7.8/I8 


2.0 


3.2±0.3 


195 


18.6/18 


1.5 


2 9+0.6 




6.87/6 


1.7 


3 1+^^-^ 

•'■^-0.5 




6.36/6 


1.8 


a O+0.6 

•'•^-o.s 




6.16/6 


2.0 


3.4 ±0.6 




5.86/6 



Note. — Spectral models are an absorbed power law 
for the faint state (Obsid 7248) and an absorbed power 
law plus narrow Gaussian, with line energy (6.4 keV) 
and redshift (0.015728) fixed, for the bright state. The 
photon index F was also not allowed to vary in the fits. 
Spectral models were fit over the 3 — 10 keV (2— 10 keV) 
for the bright (faint) state, respectively. Errors are 90% 
confidence levels. 



TABLE 3 
Mid-Infrared Photometry 
Parameters 



waveband 


aperture" 


zeropoint* 


fim 


correction 


Jy 


3.6 


1.213 


280.9 


4.5 


1.234 


179.7 


5.8 


1.379 


115.0 


8.0 


1.584 


64.1 



Note. — " Infrared Array 

Camera Data Handbook 3.0 (2006)^ 

http ://ssc . spitzer caltech . edti/irac/dh/ira cdatahandbook3.0.pdf] 
p. 53. * Vega relative zeropoints from 
Reach et al. (2005). 



10 



TABLE 4 
IC 4970 Nucleus Mid-Infrared 
Photometry 



waveband flux density magnitude 





mjy 




3.6 


4.85 


11.91 


4.5 


3.64 


11.74 


5.8 


6.53 


10.62 


8.0 


15.78 


9.02 



Note. — Results are for a fixed 
2f'44 radius circular aperture centered on 
the nucleus of IC 4970 (20'' 16'" 57 '.4, 
-70°45'00"1) with a concentric back- 
ground annulus whose (inner, outer) radii 
are (2f'44, 7f'32). Photometric correc- 
tions and zeropoints are listed in Table 
[3] IRAC photometric uncertainties are 
~ 5%. 



